Vortex domain walls in helical magnets 
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We show that heUcal magnets exhibit a non-trivial type of domain wall consisting of a regular 
array of vortex lines, except of a few distinguished orientations. This result follows from topological 
consideration and is independent of the microscopic models. We used simple models to calculate 
the shape and energetics of vortex walls in centrosymmetric and non-centrosymmetric crystals. 
Vortices are strongly anisotropic, deviating from the conventional Berezinskii-Kosterlitz-Thouless 
form. The width of the domain walls depend only weakly on the magnetic anisotropy, in contrast to 
ferromagnets and antiferromagnets. We show that vortex walls can be driven by external currents 
and in multi-ferroics also by electric fields. 

PACS numbers: 75.10.-b, 75.60, 75.70, 75.85 



Introduction. — The structure of domain walls (DWs) 
determines to a large extent the properties of magnetic 
materials, in particular their hardness and switching be- 
havior, it represents an essential ingredient of spintronics 
[U 12]. Common DWs are of Bloch and Neel types in 
which the magnetization rotates around a fixed axis, giv- 
ing rise to a one-dimensional magnetization profile [3| l4l . 
Two-dimensional vortex wall configurations can appear 
in restricted geometries as a result of the competition 
of stray field, exchange and anisotropy energy [1 . The 
more difficult problem of DWs in helical magnets has not 
yet been solved. 

Here we show that DWs in helical magnets are fun- 
damentally different from Bloch and Neel walls. They 
are generically characterized by a two-dimensional pat- 
tern. For almost all orientations of the DW they contain 
a regular lattice of vortex singularities. However DWs of 
few exceptional orientations, determined by symmetry, 
are free of vortices and maximally stable. Though DWs 
do not exist without anisotropy, their width and energy 
depend only weakly on the anisotropy strength. Simi- 
lar to other topological defects [SHE] , vortex DWs can be 
driven by electric currents. In multi-ferroics vortices are 
electrically charged, allowing manipulation of magnetic 
DWs by electric fields |9h11]. 

Helical magnets exhibit a screw-like periodic spin pat- 
tern intermediate between ferromagnets and antiferro- 
magnets. Examples of such structures are shown in 
Figjl] In addition to time reversal symmetry, in heli- 
cal magnets the space inversion symmetry is broken [12] , 
either spontaneously in centrosymmetric crystals, or en- 
forced by the symmetry of the crystalline lattice in non- 
centrosymmetric crystals. The magnetization m in these 
structures rotates around a fixed axis when the coordi- 
nate along a fixed direction, generally not coinciding with 
the rotation axis, changes. Further we denote the pro- 
jection of the magnetization to the rotation axis 7713, its 
rotating projection to the perpendicular plane as mx and 



assume that m = 1. The angle of rotation is (j). 




FIG. 1. Different types of helical ordering, (a) The mag- 
netization rotates in a plane perpendicular to the helical (x-) 
axis as in Tb, Dy, Ho. (b) Conical phase with non-zero 777,3- 
component of the magnetization as in Ho below 19K. (c) The 
magnetization rotates in a plane parallel to the helical axis as 
inTbMnOs. 

Centrosymmetric case. — We begin with the cen- 
trosymmetric case, since it is simpler and includes al- 
ready many features discussed in this article. Prominent 
experimental realizations are frustrated antiferromagnets 
in rare earth metals Tb, Dy, Ho [13l[14], their alloys and 
compounds RMnOs R G {Y, Tb, Dy} [15], R2Mn205, 
Rg {Tb, Bi}, as well as Ni3V208 and LiCu202 [EHIS]. 
The helical magnetic order originates in these materi- 
als from the indirect RKKY exchange which results in a 
competing nearest neighbor ferromagnetic (J > 0) and 
next nearest neighbor antiferromagnetic {J' < 0) interac- 
tion along the helical axis [l4l[T7[[T3. The corresponding 
Ginzburg-Landau Hamiltonian then reads [1] 
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where J^ = J d^r, V_l = ydy + zdz, and a is the lattice 
constant. 9 = arccos(J/4| J'|) denotes the angle between 



spins in neighboring layers. The continuum approach 
is vahd for ^ <C 1. can be diminished to zero under 
uniaxial pressure [19 . The last term in (IT]) is an inter- 
polation that fixes the spins either in-plane, ms = at 
T = {T - To)/To > 0, as in Tb, Dy, Ho and TbMnOs, 
or on a cone with angle i9o for r < 0, as in Ho below 
To = 19K [20] (^0 ~ 1-56 [M]). ja ^ 0.625 for Ho and 
J a = 0.17 for Tb[T4l. The ground state of (1) has a 
helical structure with (j) = qx: 

m = \m^ I (ei cos qx + x^2 sin qx) + C^m^e^ (2) 

where q = 0/a (see Figjl^). x = il and (" = ±1 describe 
the chirality and conicity of the solution, respectively. 
The rotation axis 63 may be parallel to the helical axis 
X, as in Tb, Dy, Ho, or perpendicular to it, as in TbMnOa 
(see Figjlfc). Because of its space inversion symmetry, (1) 
is a generic model for any centrosymmetric helical mag- 
net. In centrosymmetric helical magnets where the star 
of modulation vectors includes 3 vectors, like in CuCr02 
[TTJ [21], a slightly more complicated model has to be 
used, but the main conclusions of our analysis remain 
valid also in this case. 

Domain walls and vortices. — DWs separate half spaces 
with different values of (" or x or both. We consider here 
only walls with different x since domain walls between 
phases with different (", but the same value of x, are of 
Ising type and well studied. A wall whose normal fi is 
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FIG. 2. DWs in centrosymmetric helical magnets. Cross sec- 
tion parallel to the x^/— plane of (a) a Hubert wall, (b) a vortex 
wall parallel to the helical axis in a system where the magne- 
tization rotates in the x-y plane, (c) a vortex wall tilted with 
respect to the helical axis. The arrows denote the orientation 
of m. For systems where m is confined to the ^2;-plane, m 
have been rotated by 7r/2 for better visibility. The red contour 
is described in the text 

parallel to the helical axis, n • x = 1, has been stud- 
ied by Hubert [TJ |22] . In such a wall the derivative of 
the rotation phase dx(l) changes smoothly from —q to q 
over a distance ^ l/q (see Figj2k). Its surface tension 
an ^ {J/a'^)\0\^ is small for small 0. Walls of different 
orientation were not yet studied theoretically, although 
seen in experiment, e.g. in Ho by circular polarized x- 
rays [l^ . We consider first a wall in the xz-plane whose 
normal n is perpendicular to x. Since both domains have 
the same pitch, the magnetization is periodic along x-axis 
with the period 27r/g. Circulating counterclockwise along 



a closed contour C in the x?/-plane formed by two hori- 
zontal lines at X = Nn/q and x = {N -\- Ny)7r/q with N 
and Ny being integers and two vertical lines connecting 
the horizontal ones far from the wall (see the red contour 
in FigJ2]3), an observer sees the change of phase 27TNy. A 
similar contour C enclosing a Hubert wall gives A^^ = 0. 
We note that this argument is purely topological and not 
limited to the particular Hamiltonian (fTl). In the case of 
six modulation vectors ±qi,i = 1,2,3, as in CuCr02, 
in addition to the ±q^ DWs considered here, also DWs 
between qi,qj phases {i 7^ j) appear, similar to those 
discussed below for the non-centrosymmetric case. 

Vortices are saddle point configurations of the Hamil- 
tonian ([1]). For 7a ::^ 1 they obey the equation 

{4Vi + a^ [6(5,0)2 _ 2^2 _ ^2] ^2| ^ ^ Q (3) 

Vortex lines parallel to x have the standard Kosterlitz- 
Thouless form [23 . The same applies to vortex lines 
perpendicular to x on scales much larger than q~^ where 
{dx<p)'^ ~ q^ and hence eq. ([3| becomes Laplace's equa- 
tion. On smaller scales, instead of solving (|3| exactly, we 
use a variational Ansatz 0(r) = arctan(Az/x), where A 
is a variational parameter to be found from the energy 
minimization. It gives A^(r) = ^^ + 5/(64 ln(r /a)) where 



X + A 2: . The vortex energy per unit length is 



£^(r) = ^ln^/'(r/a) 
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Q describes the crossover from the conventional 
Kosterlitz-Thouless behavior ~ ln(r/a) at distances r > 
Vc = aexp[5/(64^^)] to a [ln(r/a)]^/^ behavior at scales 
r < Tc. 

So far we assumed that 7a ^ 1 and hence the 
spins are confined at a fixed value of 7713. However 
for 7a < 1 in the vortex center, i.e. for r < r^ = 

7~-^(l + rcos^'^o)""^ |1^(7<^)I 7 spins align parallel to 
the 63 -axis to save energy. Thus ms( = ±1, i.e. the 
vortex forms a meron [24]. Vortices in the DW have the 
same vorticity ±1 and are equidistant with the spacing 
ir/q forming a vortex fence. The energy per unit area of 
the vortex DW is a^ = {V5J/4a^)\0\\\n\0\\^/^ > an- 

A DW of general orientation with fi • x = cos a consists 
of a periodic chain of vortices perpendicular to the heli- 
cal axis and the normal to the DW (FigJ2fc). For a close 
to the wall can be treated as pieces of Hubert walls 
separated by vortex steps of the height n/q and length 
(tt/^)/! tana|, giving rise to a vortex staircase. The en- 
ergy per unit area of such a wall is approximately equal 
to to ey{q~^q\ sin a\/7r. At any a ^ 0, it is larger than 
the energy of the Hubert wall. 

N on- centro symmetries ease. — In these systems invari- 
ants violating the space but not time inversion symmetry 
are permitted. Those terms appear in first order per- 
turbation theory in the spin-orbit coupling constant g 
\ . Experimental examples of non-centrosymmetric 



compounds are MnSi [27^, Fei_^Co^Si ^ and FeGe [29]. 
The magnetic anisotropy in crystals with cubic symme- 
try is of the order g^. The phenomenological Ginzburg- 
Landau functional for the magnetization m has been de- 
rived in detail in 130 and takes the form 



^n = - /J(Vm)2 + 2^m(Vxm)+^^mH. 



(5) 



Here we ignored other terms representing the cubic 
anisotropy since they do not influence our results quali- 
tatively. For V = the minimum of energy ^ is given by 
a planar chiral structure, m(r) = ei cosqr + e2sinqr, 
where q is the wave vector of the helix and ei, e2 = qxei 
and q form a triad. The direction of q is arbitrary, but its 
length |q| = ^ is fixed. Contrary to the centrosymmet- 




FIG. 3. DWs in non-centrosymmetric helical magnets. De- 
tail of Figure Ig of Ref. ^29^ (center) showing two types of 
DWs in the ferromagnet FeGe, the left one includes vortices, 
the right one is vortex free. The panels are theoretically cal- 
culated DWs, right without vortices, left with vortices. 

ric helices, states with wave vectors q and — q describe 
the same magnetization reducing the degeneracy space to 
SO{3)/Z2 \31 . Cubic anisotropy pins the helix direction 
q either along one of the cube diagonals or along one of 
the four- fold axis, depending on the sign of v. DWs sepa- 
rate half spaces with different values of q. Since |^'| <C ^^, 
one could expect, in analogy with ferromagnets, that the 
DW locally represents a helical structure whose wave vec- 
tor slowly rotates pertaining its length constant. We will 
prove that such a configuration does not exist. Indeed, 
the generalization of the equation for the magnetization 
in a structure with slowly varying q is 



m(r) = ei cos 0(r) + 62 sin ^(r). 



(6) 



where ^(r) is an arbitrary function of coordinates. 
ei,e2,V^ form a right triad. The requirement of the 
constancy of the pitch implies (V^) = g^, which is 
the Hamilton-Jacobi equation for a free particle with the 
boundary conditions \/(j) ^ qi^2 at x ^ ^00. Since a 
free particle conserves its momentum, the latter cannot 
be different in two different asymptotic regions. Thus it 
is impossible to construct a DW between two different 
asymptotic values of the wave vector without changing 
its modulus between. The DW solution has a width de- 
termined by the only existing scale 1/q and the surface 
energy is independent of anisotropy v. 



DWs whose plane is a bisector of the asymptotic wave 
vectors qi and q2 do not contain vortices. They are 
analogs of the Hubert DWs. Their surface tension has 
the order of magnitude a ~ Jg/a. DWs of any differ- 
ent orientation contain a chain of vortex lines for the 
same reason as in the centrosymmetric case (see FigJ3J 
right panel). The vortex lines are located in the plane 
of the DW perpendicular to the projection of either the 
vector qi — q2 = 2q_ or qi + q2 = 2q+ onto the do- 
main plane depending on what configuration has lower 
energy. The vortex line spacings in the chain are equal 
to £± = 27r/ |n x q±|. Pictures of both vortex- free and 
vortex DWs based on a variational numerical calculations 
are shown in Fig. |3J together with the experimental fig- 
ure of FeGe [29l displaying these structures. For numeri- 
cal calculations we used ^ and the following ansatz (we 
write the answer for the first choice of sign) : 



^(r) = rq+ + nq_i^lncosh 



nr 



arctan 



tanh '02 



(7) 



where V^i = [r — n(nr)] q_ and ip2 = |n x q_| nr. The 
last term in ([t]) is the contribution of the vortex array. 
It has the asymptotics i^^i. The second term does not 
have any singularity. It corresponds to the vortex-free 
DW when n is parallel to qi — q2, i.e. when the DW 
plane is the bisector of the vectors qi and q2 . Its asymp- 
totics are ± (nr) (nq_). The asymptotic of the sum of the 
second and third terms is ±rq_. Together with the first 
term they tend asymptotically to qir above the domain 
wall and to q2r below. The only variational parameter is 
w. The surface tension of a vortex DW differs from that 
of the vortex-free DW by a factor sin/31n(l/ga), where 
P is the angle between n and q^. Apart from a nar- 
row interval of small /3, this factor is larger than one. 
Because of their higher surface tension, DWs carrying 
vortices may be unstable with respect to formation of 
a zig-zag structure formed by vortex-free DWs. Zig-zag 
structures observed in experiments with Feo.5Coo.5Si [28 
can be tentatively interpreted as arising from this insta- 
bility. The zig-zag structure is impossible in the helical 
magnets with uniaxial anisotropy since only one orienta- 
tion of the vortex free DWs is allowed. This fact together 
with low stability of vortex-carrying DWs can serve as ex- 
planation of a disordered domain structure observed in 
Ho [20 . 

DW roughening. — Roughening of DWs occurs by for- 
mation of terraces which condense at the roughening 
transition temperature [32]. For Hubert walls terraces 
are encircled by vortex rings of some length L. Since 
their energy and entropy scale as £y{L){L/a) and L/a, 
respectively, Hubert walls remain asymptotically fiat at 
increasing temperatures, slowing down their propagation. 
On the contrary, vortex walls are always rough, as seen 
also experimentally [2Ql . 

Driven domain walls. — We assume that the spin of 
a conduction electron follows adiabatically the magne- 



tization ni(r). This approximation is valid provided 
\kp — kp\ ^ q. Here kp is the Fermi momentum of the 
electrons with spin parallel or anti-parallel to m. Thus, 
electrons experience a change of angular momentum. In- 
versely, the electron current j creates a reaction torque 
on m driving the magnetic texture with a force [5HH] 

Fa = ^j/3 / {m • (Sam X dfsm) + Psfd^sm • ^^01} (8) 



The first term is the spin transfer torque [5^ 6 related to 
the Berry's curvature K^ = ec^/^^m (9/3 m x d^m) . For 
a single vortex its only non-zero component is parallel 
to the vortex lines and is given by 2t: m^C,. A weak field 
along the axis of rotation will order (" of different merons. 
The force per unit area of the DW exerted by a current 
of density j parallel to the wall due to the spin torque is 
of the order m^^C^^j /l{)^Am~'^)Nm~'^ . The second term 
results from the spin relaxation and is orthogonal to the 
first one. psf is a dimensionless coefficient which depends 
on the specific relaxation mechanism [7, 8 . The pinning 
force density due to non-magnetic impurities of density n^ 
can be estimated from the theory of collective pinning as 
JOrii/Q ^ 0{Tc/20K){ni/10^^cm-^)Nm-^, which gives 



a critical current jc ^ 6 10 Am for 



lO^^cm-^. 



Multiferroics. — In multiferroics the magnetization can 
induce the electric polarization [9] 



P = K [m(V m) — (m V)m] , 



(9) 



where Hi is some material constant. P is only non-zero 
if mx 7^ (as in TbMnOs). The vortex structure in a 
helical DW induces a ferroelectric DW, in agreement with 
experiments [33 . Hubert walls are uncharged whereas 
vortex lines carry an electric charge p = 27Thi [es x x] fi 
per unit length. This allows to move magnetic DWs by 
an external electric field. 

To conclude, we have shown that DWs both in cen- 
trosymmetric and non-centrosymmetric helical magnets 
consist of a regular array of vortex lines for almost all 
orientations except of a few that correspond to a minima 
of the surface energy. The helical DWs are generically 2- 
dimensional textures. They are charged in multi-ferroics 
and can be driven by electrical currents and fields. 
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